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Abstract 

The fundamental problem of Farley-Buneman turbulence in the auroral £?-region has 
been discussed and debated extensively in the past two decades. In the present paper 
we intend to clarify the different steps that the auroral E'-region plasma has to undergo 
before reaching a steady state. The mode-coupling calculation, for Farley-Buneman 
turbulence, is developed in order to place it in perspective and to estimate its magni- 
tude relative to the anomalous effects which arise through the nonlinear wave-particle 
interaction. This nonlinear effect, known as nonlinear "Landau damping" is due to 
the coupling of waves which produces other waves which in turn lose energy to the 
bulk of the particles by Landau damping. This leads to a decay of the wave energy 
and consequently a heating of the plasma. An equation governing the evolution of the 
field spectrum is derived and a physical interpration for each of its terms is provided. 

1 Introduction 

The selfconsistent theory of Farley-Buneman turbulence as developed by Sudan [8, 9], 
and by Hamza and St-Maurice [3, 4] has addressed successfully a number of features 
observed by radars in the equatorial and auroral ionosphere. Recently, Hamza and St- 
Maurice [5, 6] developed a general fluid formalism, taking into account the previously 
ignored nonlinear electron inertia effect term. The new turbulence theory provides a 
possible explanation for the existence of large aspect angle echoes in the auroral E'-region, 
and justifies selfconsistently the link between electron heating along the magnetic field 
lines and the anomalous effects induced by perpendicular Farley-Buneman turbulence. In 
[5, 6] it is shown that the advection of parallel current by perpendicular drifts leads to 
an anomalous parallel collision frequency. In other words, it is shown that perpendicular 
Farley-Buneman turbulence can indeed lead to anomalous parallel effects which in turn 
allow for parallel electric fields to be generated. This has an impact on energy conservation 
since the zeroth order Ohmic heating which was confined to the plane perpendicular to 
the magnetic field, can now have a parallel component due to plasma wave turbulence 
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(J e ||.E||). This is made possible by the nonlinear electron inertia term. In order to clarify 
the physics related to the energy transport properties of the £ , -region plasma we have to 
discuss the different mechanisms of energy tranfer from the free energy source confined 
to the plane perpendicular to the magnetic field to the particles moving along the field 
lines. This can only be achieved by including mode-coupling effects and describing their 
role as it pertains to the final picture. 

In what follows, we will use the results of Hamza and St-Maurice [5, 6] to establish 
the final link by identifying the mode-coupling terms and the terms associated with the 
analog of a wave-particle interaction respectively. In particular, we will stress the fact that 
the mode-coupling terms are energy conserving and that they do contribute a frequency 
broadening Auj^ as already shown by Hamza and St-Maurice [3] as well as Sudan [8, 9]. 
Finally we will discuss the roles of the mode-coupling terms and that of the anomalous 
effects generated via nonlinear wave-particle analogs respectively. However, more impor- 
tant than the comparison of the relative magnitudes of the mode-coupling rate and the 
parallel rate of dissipation generated by the anomalous effects is the difference in nature 
of the physics associated with each respectively. Although, the frequency broadening gen- 
erated by mode-coupling, as shown by Hamza and St-Maurice [3], might indeed balance 
the linear growth rate of a Farley-Buneman wave, it is only shuffling energy into other 
waves. The parallel dissipation rate, on the other hand, allows wave energy to be chan- 
nelled directly into electron kinetic energy. In steady state turbulence a physical process 
must ultimately convert wave energy into thermal energy. This energy conservation, in 
fact is a requirement for any selfconsistent theoretical model. 

We already know that during the so-called "quasi-linear" phase, the background particle 
distribution diffuses in such a way as to bring the growth rate of the unstable waves to zero 
( "plateauing" of the distribution function), leaving a quasi-stationary spectrum interact 
in such a way as to distort the spectrum but keep the energy in the spectrum roughly 
constant. This is achieved by extracting energy from the resonant particles and dumping 
it unto the waves. Quasi-linear theory leads to the development of a quasi-equilibrium 
spectrum which persists indefinitely. The spectral energy can be shown to be proportional 
to j/uj, where 7 and u represent the growth rate and the eigenfrequency respectively. This 
energy is small since it is proportional to j/u <C 1, and therefore can not be considered 
to be an ultimate saturation mechanism for the problem at hand as shown by St-Maurice 
[7]. Mode-coupling effects, neglected by the quasi-linear theory lead to a broadening of 
the equlibrium spectrum. The mode-coupling effects can be divided into two categories. 
The first class of terms can be qualified as "resonant" in the sense that waves couple while 
satisfying the following relations (frequency and wavemomentum matching) 

w k'+ w k" =a; k, 

k'+k"=k. [ ' 

The resonant mode coupling terms lead to a broadening of the quasi-linear spectrum and 
the time scale associated with this effect is proportional to the spectral density (in our 
case \8(f)\ 2 oc 7), and therefore the resonant mode coupling process is on the same time 
scale as the quasi-linear time scale. However, as shown by Hamza and St-Maurice [3, 4] 
and Drummond and Pines [2] before them, the resonant mode coupling leaves the total 
energy in the wave spectrum virtually untouched (constant). 



464 



A.M. Hamza 



One major constraint in resonant mode coupling is expressed by equation (||). The 
driving term is a quadratic nonlinearity of the form ^'W^k-k'^) an d induces a time 
dependence of exp [i(w k ' + w k _ k ')i], which in general can be very different from the time 
dependence of (ftk(t), that is, exp [iw^t]. In other words, in the general case one has 

w k' + "V + w k, 

k'+k"=k. [) 

The nonlinear source forces the field <j>k(t) to have terms with the time dependence of 
the source rather than the natural time dependence of the field itself. This is very much 
the case for a forced harmonic oscillator. This leads to separating the field variable into 
two parts. The first part gives rise to resonant coupling because of frequency matching, 
while the second part arises because of the non-resonant parts of the nonlinear source. 
An essential feature of the latter part is that (uv/ + w k ») can be fairly small so as to 
allow the phase velocity of the resultant wave (w k ' + ^ k _ k ')/^|| to be smaller than the 
electron thermal velocity, and for that matter the ion thermal velocity. This means that 
the non-resonant contribution of the mode-coupling terms can undergo Landau damping. 
The non-resonant waves Landau damp and in the process give their energy to both the 
electron and ions, i.e. heat the plasma along the field lines. 

The main objective of the present calculation, is to show that the nonlinearities in 
Farley-Buneman turbulence contribute coherent and resonant terms. These terms are 
isolated and shown to have different impacts on the physics, especially on the energy 
conservation constraint. We have omitted to discuss transient phenomena because we 
are only interested in the global energy conservation on a long time scale. Higher order 
nonlinearities have also been neglected in the perturbation analysis. The model has been 
extensively described in [6]. Therefore we will refer the reader to that paper for formal 
manipulations and for details related to the derivation of certain fundamental results. 



2 Theoretical Model 

The fundamental assumptions and the details of the theoretical model are described, as 
mentioned above in [6]. We will only quote the results neccessary for the development of 
an energy balance equation containing all the relevant terms including the mode-coupling 
terms. 

The starting point of this paper is equation (24) of [6] which can be written in a compact 
form as follows 

-zn ^(^-4 L) )=( 1 ) + ( 2 ) + ( 3 )' (3) 
where the three terms on the right hand side of equation (|3]) can be explicitly written as 

{) - l + Vwe^ kj A k ,fAk' T e J ek-kW W 

k w 

This term represents the advection of parallel current by perpendicular and parallel drifts 
through electron inertia. The second term corresponds to the advection of perpendicular 
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current by parallel and perpendicular drifts and can be written in the following form 
lJ -l + Vk*e 2 k\ ^ + AL 

(5) 

x V ( Qe - ^ I A^ J \icrv ■ J- 1 - , , 

k'w' v k-k'w-o/ " e / Je 

It is important to emphasize the fact that these two terms were completely omitted by 
previous theories addressing the issues related to Farley-Buneman turbulence in the auroral 
as well as equatorial E region. 

Finally the last term represents the classical mode coupling terms, which can be written 
in the following form 

< 3 > - - £ w H» - ^ K^) } 



\- W ,a,;q/ • (k x k') e0^y 

a . _ _ , V C k;k' |k , - k' , I T P ^-k'a,-a,' • 



i ■' ' k— k a;— a; 

K 

In the above expressions we have 

while a> k represents the linear eigenfrequency for Farley-Buneman waves. The coefficients 
A'fl'f,, ti?' 1 ?,, and C?] 1 ?, are defined in the paper [6], and the parallel and perpendicular 

currents jj and are given by equations (21) and (22) in [6] respectively. 

It is clear from equation (Q) that when the nonlinear right hand side is neglected we 
recover the classic Farley-Buneman dispersion relation. 

The first two terms on the right hand side of equation (]3|) contribute the coherent terms 
discussed extensively in the paper [6]. The coherent terms can easily be identified through 
a formal substitution of the expressions for jj and J^r into equation @. Isolating these 
coherent terms, while neglecting the four- wave coupling which generates cubic nonlinear- 
ities, allows us to rewrite the top equation @ in the following form. To illustrate this 
iteration process let us consider the first term in equation (|3|). To proceed we need the 
expression for the parallel current 



|| .n e 2 fc|| x ^ \y || 



(8) 



= 2 „ k n 
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Substituting equation @ into equation (|4|) one obtains two types of terms. The first type 
can be combined with the mode coupling term given by equation (|B|), while the second 

type of term, a cubic term, is of the form (jxfiJe. A second iteration of this term leads 
to a cubic term of the form \(j>^j | 2 0kw This term can combined with the left hand side 
of equation (|3|) to obtain a renormalized dispersion relation. This effect introduces an 
anomalous frequency solely due to nonlinear effects analog to a higher order wave-particle 
interaction. 

- i(w - w k )(^ = M k ; k' ^jK-^u-J . ( 9 ) 

1 ' ' 

k uj 

where M^'£, is the nonlinear coupling coefficient obtained by substituting the expressions 

for jj and J^, as given in [6], in equation @ above. Where the new frequency Wk includes 
the coherent terms generated by the nonlinearities, and written as i/HjJ and v^J-^, where 
jj and 3jr are shown to be propotional to (f)^ to first order. These nonlinear coherent 
expressions were derived and discussed extensively in [6]. a>k can now be written as follows: 

w k = 4 L) + ^k , (io) 

where ui^ represents the anomalous (coherent) contribution by the nonlinear wave-particle 
analogs in equation (|J). As one can expect, the nonlinear expression for this anomalous 
frequency is quite complicated and will not be given here. 

Making the following change in wavevectors and frequencies k — > k — k , and uj — ► 
uj — uj leads to rewriting equation (||) in a useful form for our purpose as will become clear 
in what follows. 

- i{0J - CJk)^ = ~ ( M k;k' + M k;k-£ ) ^kV K-V' U -J ■ ( n ) 

k'a/ 

Multiplying by and ensemble averaging leads to the evolution equation for the spectral 
density (l^l 2 ) 

" i(« " ^k)(|^| 2 > = \Y (Kk' + K;k^') (^kV^k-kW^L) ■ ( 12 ) 

k u> 

At this stage one has to close the system by prescribing the triple correlation function. We 
shall assume that the distribution of potential fluctuations is Gaussian, which allows to 
evaluate explicitly, via a perturbation scheme, the right hand side of equation (|l2|), that 
is 

(j) = + 0( 2 ) + • • • 

with 

-i(.-4 L, )« = ». 

- = \ £ (<£ + *£X- t w <13> 

k uj 
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the triple correlation function can be written as: 

(14) 

The first term on the right hand side of equation (Il4|) vanishes, while the others can be 



combined with equation (13) to lead to 



,(1) ^(1) ^(2)*\ _ 2 ^ V k;k k;k-k j „ „ 

2 , , , - ^k) 



, +M, , , 



^(2) _ i V V k-k;-k k-k;k; (15) 



y ku 



, , / , + M , 



2 , , , (^ - w k' ) 



k' I 
u> 

This finally allows us to write the equation (12) for the spectral density in the following 



form: 



- -k)(i0 k j 2 ) = ± E + Kk-k'') 

k a; 

+ M kik ) <I<W l 2 )d^l 2 > (is) 



(w-w' -w k _ k ') 



1 

~ 2 (w - w k ) V^ k;k ' ' ^ k ; k - k ' 

The first term on the right hand side corresponds to a broadening due to resonant three- 
wave coupling (mode coupling), and can be incorporated with the left hand side to give 

-i(uj - Wk - Aw k )(|</> ka; | 2 ) 

i^rj E + M k Tk'T <i^v i 2 xi^- k w i 2 >> (17) 

k U! 

where 



M k;k' +M k;k-k' J l iW k-k';-k' +M k-k';k , 

^ k " "2 ^ 7^ V ^ ^l^kVl 2 )- (18) 

k w 

These results are very much similar to the two dimensional results published by Hamza and 
St-Maurice [4] when looking explicitly at the two dimensional Farley-Buneman turbulence. 
The coefficients in the current case are more complicated because we have taken into 
account the electron inertia nonlinearities and included the parallel effects. We do not 
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intend to explicitly write the expressions for the coefficients since, as might be expected 
they are quite complicated and would not serve the purpose of this paper. However, we 
have in equation ( |17| ) isolated two types of terms. The first term is a coherent analog of 
a nonlinear wave-particle interaction, nonlinear Landau-damping, while the second is due 
to a resonant three-wave coupling. The term on the right hand side of equation ([I?]) is 
due solely to mode coupling. 

At this stage, one can study two limiting cases, namely a strong turbulence limit and 
a weak tubulence one. 

2.1 Strong Turbulence Limit 

In this case equation ( |l7|) can be written in a more familiar form: 

,x * 



k 10 

This equation can be written in the following form 

(|^| 2 ) = ^^, (20) 

where (|0k^| 2 ) is the incoherent driven potential obtained by identifying the right hand 
side of equation (|20|) to the right hand side of equation (|l9|). This has been the object of 
a long debated issue related to describing strong turbulence as it occurs in plasmas. We 
will try to illustrate, very briefly, the results of the debate on strong turbulence using the 
fundamental results of Boutros-Ghali and Dupree [1] and references therein. When writing 
Poisson's equation we separate the coherent (phase coherent with the electric field) from 
the incoherent contribution of the distribution function in the following way 

V • 5E = Annq J dv (5f c + sfj , (21) 

where Sf c and 8f represent the coherent and incoherent parts of the particle distribution 
function respectively. The coherent part contributes to the dispersion relation which allows 
us to write 

ie kuj k • 5E = Airnq J dvdf (22) 

using <5E = — ik</>ku leads to 

47rng f , ^ 

3kw = — / dvdf = (23) 



which basically leads to equation ( p0| ) . A strong turbulence result (see for example Boutros- 
Ghali and Dupree [1]). 
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2.2 Weak Turbulence Case 

In the weak turbulence regime, in which the waves are weakly growing, with small line 
widths Im (w^), and replacing Fourier series with Fourier integrals, we obtain 

£<l*U a > - / (^4 / (24) 



and using the following spectral relation 



2 'f) k 27rd(u-u k ). (25) 



/kcj 



Let us now go back to equation (17) and integrate over the frequency using the following 
identity 

- 2Re J ^i(co - 4 L) XI0U 2 ) = |(<A 2 )k. (26) 

In other words, the time derivative represents the quasi-linear growth rate. We finally 
obtain 

d , 2 . - _ f du> 



2Re J ^(u£ + Aw k )(|0 k J 2 ) 

f dk' / (27) 

~ 2R(i J (2^)3 + M k;k-k'J (^)k' <</>Vk'^k - U k , - U k _ k ,) 

the first two terms contribute two different growth rates: 7^ due to the coherent nonlinear 
wave-particle analog, and 7^ due to mode coupling effects. This allows us to write the 
following spectral equation 

|(0 2 ) k = -2(^ + ^)(0 2 > k 

f dk' , - (28) 

~ 2Re J ( Mk *' + M k;k-k'J <^>k' <0\-k'^k - <*V - u, k _ k ,) 

integrating over k one can show that the mode-coupling term conserves the total mode 
energy, which leaves us with: 

dk|(0 2 ) k = -2 J dk( 7k 7 )(0 2 ) k (29) 

in other words, in order to reach a steady state the right hand side has to vanish. This 
leads us to the following interpretation. A steady state requires both resonant mode 
coupling and the analog of a wave-particle exchange to take place; one in the absence of 
the other does not neccessarily guarantee steady state turbulence. Consequently, a steady 
state in Farley-Buneman turbulence requires that the free energy available in the plane 
perpendicular to the magnetic field, extracted via a linear instability, be converted into 
particle kinetic energy along the field lines via an analog of a nonlinear wave-particle 
mechanism. 
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3 Summary and Conclusions 

In this letter we have established the neccessary condition for steady state Farley-Buneman 
turbulence. We have shown that in order to reach steady state the growth in the plane 
perpendicular to the magnetic field has to be compensated by a growth along the field 
line; A perfect analogy with a leaky bucket, if the we want the water level in the bucket to 
remain unchanged then we need to balance the rate of pouring to the rate of leaking water. 
In other words free energy is extracted from the drifts by the Farley-Buneman waves which 
start to grow and couple. This process will go on, and no steady state will be reached 
unless the energy is given back to the particles, in this case the particles moving along 
the field lines. So two processes take place. First, linear growth and resonant coupling 
of Farley-Buneman waves, then a nonlinear wave particle coupling along the field lines. 
These two mechanisms when combined together can lead to steady state turbulence. 
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